High accuracy calculation of 6s — ► 7s parity nonconserving amplitude in Cs. 



M. G. Kozlovf] and S. G. Porsev 
Petersburg Nuclear Physics Institute, Gatchma, 188300, Russia 



o 
o 

(N 

Oh- 
< 

in 

(N 



Oh: 



c3 

C/3 

o 

C/3 . 

Oh: 



> 
m 
in 
o 



o : 

(Z3 1 

>>: 
43 : 

Oh' 



X 



I. I. Tupitsyn 

Institute of Physics, St. Petersburg State University, 
Uljanovskaya 1, Petrodvorets, 198504, St. Petersburg, Russia 
(Dated: February 2, 2008) 

We calculated the parity nonconserving (PNC) 6s — > 7s amplitude in Cs. In the Dirac-Coulomb 
approximation our result is in a good agreement with other calculations. Breit corrections to the 
PNC amplitude and to the Stark-induced amplitude f3 are found to be —0.4% and —1% respectively. 
The weak charge of 133 Cs is Qw = —72.5 ± 0.7 in agreement with the standard model. 
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Introduction The last and the most accurate mea- 
surement of the ratio of the parity nonconserving (PNC) 
amplitude in Cs to the Stark-induced amplitude (3 was 
done in Boulder |5J]. Later the amplitude (3 was mea- 
sured with a very high accuracy 0. These measure- 
ments together with atomic calculations of the PNC am- 
plitude H Q gave the most accurate atomic test of the 
standard model. The authors of the the calculations |3|, ^] 
estimated their accuracy to be 1%. After that Bennet 
and Wieman 0] compared theoretical and experimental 
results for a number of observables including hyperfine 
constants, El amplitudes, and the Stark-induced ampli- 
tude [3 and suggested that the actual accuracy of atomic 
calculations of the PNC amplitude was 0.4%. Their anal- 
ysis gave the following value of the weak charge of 133 Cs: 
Qw = —72.06(44), which differs by 2.3a from the stan- 
dard model prediction [H: 



-73.09(3) 



This result caused an intensive discussion (see Pj and 
references therein). It also stimulated a new wave of 
atomic calculations 0, ||, First Derevianko H sug- 
gested that previous calculations underestimated Breit 
correction. He calculated the latter and added it to the 
PNC amplitude from 13, 0]. That reduced the discrep- 
ancy to l.Ocr. Somewhat smaller Breit corrections were 
obtained in ||, |j| . 

Here we re-analyze the role of Breit interaction in cal- 
culation of the PNC amplitude in Cs. In contrast to the 
paper we made a full scale calculation which (partly) 
accounts for the higher order corrections of the many- 
body perturbation theory (MBPT). The Breit interac- 
tion was included on all stages of the calculation, i.e. 
on the stage of the solution of the Dirac-Fock equations, 
then on the stage of the solution of the random phase ap- 
proximation (RPA) equations, and, finally, on the stage 
of the evaluation of the self-energy operator. All these 
corrections are of the same order of magnitude, which is 
not surprising taking into account that both RPA and 
MBPT contributions to the answer are quite significant. 



Our final value of Breit correction to the PNC amplitude 
is —0.4% in agreement with our preliminary result ||. 

We also calculated Breit correction to the Stark- 
induced amplitude /?6s,7s and found it to be about 1%. 
That suggests that an excellent agreement between the 
theory and the experiment Q for this constant may 
be partly accidental. In the end of this paper we discuss 
corrections, which are not included in PNC calculations 
and conclude that the accuracy of the atomic theory is 
close to 1%. Within this accuracy there is no contradic- 
tion with the standard model. 

Details of the calculation The Dirac-Fock equations 
were solved on a radial grid with the help of the code Jicf] 
and all MBPT calculations were done using summations 
over the basis set. We formed the Dirac-Fock-Coulomb 
(DFC) and the Dirac-Fock-Coulomb-Breit (DFCB) vari- 
ants of the basis set and repeated all calculations in ex- 
actly the same way for both of them. 

It is known from previous calculations of PNC effects in 
Cs that high orders in the residual Coulomb interaction 
are very important. Indeed, the MBPT effects account 
for 13% of the ionization potential, for more than 60% 
of the hyperfine constants, and for 22% of the PNC am- 
plitude. Here we calculated explicitly only the second 
order diagrams but introduced screening coefficients Ck 
for the Coulomb lines in these diagrams jy]. The latter 
can be calculated as an average screening of the two elec- 
tron Coulomb radial integrals of a given multipolarity k. 
These coefficients served as an approximation to the in- 
sertion of polarization operator in Coulomb lines. The 
coefficients Ck were chosen as follows: Cq = C\ = 0.85, 
and Ck = 1 for k > 2. With these screening coeffi- 
cients the Brueckner energies of the orbitals appeared 
to be within 1% from the experimental values. Still the 
additional corrections were necessary in order to reach 
the accuracy better than 1%. Thus, we used the energy 
dependence of the self-energy operator a(e) to improve 
the agreement between calculated energy and the exper- 
iment. The optimal one-electron energies e appeared to 
be close to the Brueckner energies of corresponding or- 
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bitals. In a similar way we chose the energies in the RPA 
equations for hyperfine and PNC interactions. 

The optimal choice of the onc-clcctron energies for 
<j(e) and for RPA matrix elements allowed us to improve 
the accuracy for the lower part of the spectrum and for 
the hyperfine constants. Corresponding corrections were 
smaller than 1% for the energies and 1% - 2% for the 
hyperfine constants. RPA corrections to the El ampli- 
tudes were much smaller and we did not do any fitting for 
them. The RPA equations for weak interaction Hp were 
solved for the same energies as were used for hyperfine 
constants. That increased PNC amplitude by 1%. 

All observables were calculated using Brueckner or- 
bitals and RPA amplitudes. We also calculated two 
smaller MBPT corrections to the amplitudes. One of 
them is the so called structural radiation and the other 
is the normalization correction j^, . The former is typ- 
ically of the order of few percent and the latter reduces 
the final answer by approximately 1%. There are sev- 
eral corrections to the PNC amplitude of the same order 
of MBPT as these ones. They correspond to the RPA- 
type diagrams and to the structural radiation-type dia- 
grams with both El and PNC vertexes on the inner lines. 
These corrections were shown to be negligible || || and 
we omitted them here. 

Breit interaction for valence electrons The frequency 
independent Breit interaction between electrons 1 and 2 
has the form 



V B = V£ 
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(1) 



where ati are Dirac matrices and r\i is the distance 
between the electrons. The operator Vg is called the 
Gaunt term and the operator Vg is called the retardation 
term. For valence electrons the Gaunt term is known to 
be about an order of magnitude larger then the second 
one |l2| . The latter is of the same order as the frequency 
dependent corrections to the operator (|l|) and as other 
QED corrections [12]. Therefore, below we neglect the 



second term in Eq. ( 

Calculations of Cs are done in the V N_1 approxima- 
tion and one needs to calculate Breit interaction between 



the electron and the closed shells only. It is easy to show 
that in this case only the exchange term of the Breit in- 
teraction does not turn to zero. This is almost obvious 
because for a closed shell (a) = 0. The exchange Breit 
interaction can be included in calculation cither pertur- 
batively or self-consistently. It is known jl2[ [l4| that the 
core relaxation results in a screening of Breit interaction 
between the innermost core shells and the valence elec- 
tron, significantly reducing the final Breit corrections. In 
the self-consistent approach, the potential Vg ore should 
be included in the Dirac- Fock equations. For this pur- 
pose we modified the Dirac-Fock code [|l0| to allow for 
the Coulomb-Breit potential. 

The self-consistent approach for such observables as 
hyperfine constants or PNC amplitudes requires the so- 
lution of the RPA equations. There are two types of 
Breit corrections to these equations. First, one has to 
use DFCB orbitals and DFCB orbital energies in RPA 
equations. Second, it is necessary to include Breit inter- 
action explicitly in two-electron matrix elements. These 
types of corrections are of the same order of magnitude 
and should be included on equal footing. 

Finally, we calculated the the self-energy operator also 
using DFCB orbitals and orbital energies. Here we did 
not include Breit correction to two-electron matrix ele- 
ments in diagrams for the self-energy. These diagrams 
arise in the second order in residual Coulomb interaction 
and always include two Coulomb lines. If in one of these 
lines the Coulomb interaction is changed to the Breit 
one, the resultant correction appears to be very small. 
The reason is that main Coulomb corrections are asso- 
ciated with the interaction of the valence electrons with 
the upper core shells, while the most important Breit 
interaction is with the innermost core shells. 

All calculations with the exception for structural ra- 
diation and normalization corrections were repeated two 
times in exactly the same manner in Dirac-Coulomb and 
Dirac-Coulomb-Breit approximations. That allowed us 
to single out Breit corrections even when they were com- 
parable to the errors caused by the incompleteness of the 
basis set or by the approximations we used to account 
for higher order terms of MBPT. 

Numerical results The PNC amplitude of the 6s — > 7s 
transition can be written as 



i-i ia n \ / (7sj,m\El\npj,m)(npj,m\Hp\6sj 1 m) (7s j ,m\H P \np j ,m)(np j ,m\El\6s j ,u, / , 
i/lpNclos, 7s) — \ ; ; 1 ), (2) 



n 



where j = 1/2. This amplitude is sensitive to the energy 
spectrum and .El-amplitudes as well as to the matrix el- 
ements of the weak interaction Hp. The latter depends 



on the wave function in the vicinity of the nucleus and 
in this respect is similar to the matrix elements of the 
hyperfine interaction. Thus, in order to estimate the ac- 
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TABLE I: Binding energies for the lower levels of Cs (au). 
The columns DF and BA correspond to the Dirac-Fock and 
Brueckner approximations. The column OE corresponds to 
Brueckner approximation with self-energy calculated at 'op- 
timal' energies e(6s) = —0.22, e(7s) = —0.15, and e(np) = 
—0.10. Finally we give Breit corrections to the energy in 
Brueckner approximation. 



TABLE III: PNC amplitude _El PN c(6s, 7s) in the units 



DF 



BA 



OE 



Breit Exper. 



6si/ 2 

6pi/2 
6f>3/2 
7Sl/2 

7^3/2 



0.12737 
0.08562 
0.08379 
0.05519 
0.04202 
0.04137 



0.14424 
0.09226 
0.08967 
0.05869 
0.04389 
0.04305 



0.14312 
0.09217 
0.08959 
0.05864 
0.04386 
0.04303 



0.00002 
-0.00004 
-0.00001 

0.00000 
-0.00001 

0.00000 



0.14310 
0.09217 
0.08965 
0.05865 
0.04393 
0.04310 



TABLE II: Hyperfine constants for 133 Cs (MHz). 



DF 



MBPT 



OE Breit Theory Exper. 



6«l/2 

6pi/2 

6P3/2 
7*1/2 
7P1/2 
7P3/2 



1424 
160.9 
23.9 
391.4 
57.6 
8.6 



2350 

298.6 
53.4 

546.8 
94.2 
17.1 



-52 
-5.3 
-2.3 
-0.8 
-0.2 
-0.0 



+5.0 
-0.2 
-0.0 
+0.8 
-0.0 
-0.0 



2302 

293.5 
51.2 

546.8 
94.0 
17.1 



2298 
292 

50.3 
546 

94.3 



curacy of the calculation of the PNC amplitude one has 
to analyze the accuracy for all these observables. 

Tables |j] and |l| present our results for the energies and 
hyperfine constants of the lower levels of Cs. These ta- 
bles show results on the different stages of calculation in- 
cluding the initial Dirac-Fock approximation, the MBPT 
and Breit corrections. We also give corrections associ- 
ated with the optimization of the energies in MBPT and 
RPA calculations as described above. 

Breit corrections to El amplitudes appear to be mostly 
smaller than 0.1%. Similar corrections to scalar polariz- 
abilities of the levels 6s and 7s are about —0.1%. These 
polarizabilities can be used to calculate the Stark shift 
8vq s j s for the 6s — > 7s transition (in Hz/(V/cm) 2 ): 



Su 6 



s,7s 



.7268 without Breit, 

.7259 with Breit, 

.7223 theory [4], 

.7262 experiment lljJ 



(3) 



Breit correction to the Stark-induced vector amplitude 
06s,7s is more important because this amplitude turns 
to zero in the non-relativistic approximation. Numerical 
values for this amplitude in au are: 



27.17 without Breit, 

26.89 with Breit, 

27.0(2) theory [4], 

27.02(8) experiment pi 



(4) 



In the experiment p5[ the amplitude (3 was measured rel- 
ative to the hyperfine Ml amplitude. Dzuba and Flam- 
baum suggested that the latter is slightly smaller 



10" 



)w/{-N) au. 





DF 


MBPT OE 


Breit 


Total 


This work 


-.742 


-.896 -.905 


.004 


-.901 


Dzuba et al. 










-.908 


Blundell et al. \ 






-.907 


.002 


-.905 


Derevianko 








.008 




Dzuba et al. 




] -.739 


-.907 


.005 


-.902 



(H and 



gave the value = 



than was assumed in 
26.96(5). 

Our results for the PNC amplitude are presented in ta- 
ble III . Calculations were done for the uniformly charged 
nucleus of the radius 6.20 fm. Following [Q, [l7| we intro- 
duced —0.1% correction to account for the difference in 
neutron and charge radii of the nucleus. 

We see that ab initio calculation and calculation with 
'optimal' energies for the self-energy operator and for 
RPA amplitudes differs by 1%, i.e. PNC amplitude ap- 
pears to be less sensitive to the variation of the energies 
used for the self-energy operator a(e) and for RPA equa- 
tions than the hyperfine constants. Thus, we assume that 
uncertainty for the PNC amplitude associated with the 
fitting is smaller than for the hyperfine constants. 

Discussion For the Dirac-Coulomb approximation 
our value for the PNC amplitude is in a good agreement 
with the results of the calculations |3[ Q . Our value of the 
Breit correction is approximately two times smaller than 
the result of Derevianko. This discrepancy is almost elim- 
inated if we add —0.003 correction to his result, which is 
caused by the Breit contribution to the energies of atomic 
levels. Derevianko argues that this correction should not 
be included if one fits the energies semiempirically. 

We first want to note that calculations ||, did not use 
any energy fit. Therefore, Breit correction to these cal- 
culations should definitely include this contribution. In 
our calculation we did use the energy fit and thus Dere- 
vianko's argument seems to work here. However, we did 
not insert experimental energies into Eq. (||). Instead, 
we modified the MBPT by changing the argument of the 
self-energy operator a(e). That changed the energies and 
the Brueckner orbitals thus affecting the numerators as 
well as denominators in Eq. (||). We think that this is a 
more consistent way to account for the high orders of the 
MBPT. We repeated the energy fit after Breit interaction 
was included. That changed Breit corrections to differ- 
ent contributions to the PNC amplitude, but the overall 
correction still remained —0.4%. 

The analysis of the numerical data from Tables § [] 
and Eqs. (||) and (|j) shows that the energies, the Stark- 
shift 8v, and even the Stark-induced amplitude (3 are in 
a very good agreement with the experiment. However, 
the accuracy for the hyperfine constants of the levels 6s 
and Qpi/2 is significantly lower. In order to improve these 
constants we had to make the energy fit which changed 



4 



both constants by approximately 2%. The reason for the 
low accuracy here is obvious: the MBPT corrections to 
A(6s) and A(<op 1 / 2 ) are 65% and 81% correspondingly. 
Taking into account that MBPT correction to the PNC 
amplitude is only 22% and that the energy fit changes 
this amplitude only by 1%, we estimate our accuracy for 
the PNC amplitude to be about 0.5%. 

This estimate does not include uncertainties associated 
with the nuclear structure and QED corrections. QED 
corrections to the hyperfine structure and the uncertain- 
ties associated with the nuclear magnetic moment dis- 
tribution were discussed in |ll| [ll], 20 1. The authors 
conclude that there are several corrections on the scale 
of 0.1-0.3%. 

For the PNC amplitude the vacuum polarization cor- 
rection was estimated in |2l], |2^] to be 0.2-0.4% for 
Z = 55. The core relaxation can screen it significantly, 
as it does with the Breit correction, so it is rather an 
order of magnitude estimate of the QED corrections to 
the PNC amplitude. Another uncertainty of the order of 
0.1-0.2% is caused by the poor knowledge of the neutron 
distribution in the nucleus |i~7j ]. 

In our analysis we use Gaunt approximation to Breit 
operator while in the paper j^] the retardation term was 
also included. This may explain the difference between 
our value 0.004 and the value 0.005 obtained in B. As 
we mentioned above, the retardation term is of the same 
order of magnitude as the frequency-dependent correc- 
tions to the Breit operator and as other QED correc- 
tions. Therefore, the difference between our results can 
serve as an estimate of the accuracy of Breit approxima- 
tion. Thus, the total uncertainty for the PNC amplitude 
is about 1% and there is no meaningful difference with 
the standard model. Our final value for the weak charge 



Qw = -72.5(7) 
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